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In  this  work,  electrical  capacitance  tomography  (ECT)  and  neural  networks  were  used  to  automatically  identify 
two-phase  flow  patterns  for  refrigerant  R-134a  flowing  in  a  horizontal  tube.  In  laboratory  experiments,  high¬ 
speed  images  were  recorded  for  human  classification  of  liquid-vapor  flow  patterns.  The  corresponding 
permittivity  data  obtained  from  tomograms  was  then  used  to  train  feedforward  neural  networks  to  recognize 
flow  patterns.  An  objective  was  to  determine  which  subsets  of  data  derived  from  tomograms  could  be  used  as 
input  data  by  a  neural  network  to  classify  nine  liquid-vapor  flow  patterns.  Another  objective  was  to  determine 
which  subsets  of  input  data  provide  high  identification  success  when  analyzed  by  a  neural  network.  Transitional 
flow  patterns  associated  with  common  horizontal  flow  patterns  were  considered.  A  unique  feature  of  the  current 
work  was  the  use  of  the  vertical  center  of  mass  coordinate  in  pattern  classification.  The  highest  classification 
success  rates  occurred  using  neural  network  input  which  included  the  probability  density  functions  (in  time) 
for  both  spatially  averaged  permittivity  and  center  of  mass  location  in  addition  to  the  four  statistical  moments 
(in  time)  for  spatially  averaged  permittivity  data.  The  combination  of  these  input  data  resulted  in  an  average 
success  rate  of  98.1%  for  nine  flow  patterns.  In  addition,  99%  of  the  experimental  runs  were  either  correctly 
classified  or  misclassified  by  only  one  flow  pattern. 

©  2015  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Understanding  the  two-phase  flow  behavior  of  refrigerants  is  im¬ 
portant  for  the  design  of  advanced  aircraft  cooling  systems.  Specifically, 
the  classification  of  liquid-vapor  structures  into  flow  patterns  is  useful 
for  predicting  heat  transfer  rates  and,  ultimately,  system  performance. 
Most  flow  and  heat  transfer  correlations  require  a  priori  knowledge  of 
the  two-phase  flow  pattern  and  are  based  on  steady-state  conditions 
[1,2].  Although  flow  pattern  identification  can  be  performed  using 
high-speed  imaging,  this  method  generally  relies  on  the  visual  interpre¬ 
tation  of  liquid-vapor  patterns.  Unfortunately,  visual  interpretation  can 
be  highly  subjective  [2].  As  a  consequence,  numerous  flow  pattern 
classifications  have  been  defined  in  the  past  [3-7].  In  contrast  to  visual 
observation,  which  is  often  impractical,  non-visual  sensor  signals  can 
be  analyzed  to  provide  more  objective  classifications.  Moreover,  sensor 
signals  are  desired  as  inputs  to  real-time  modeling  and  control. 

Capacitance  techniques  are  non-invasive  and  rely  on  differences  in 
electrical  permittivity  to  distinguish  between  liquid  and  vapor  phases 
[2,8-9].  Capacitance  measurements  acquired  simultaneously  with 
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high-speed  videos  of  air-water  flows  have  been  used  to  identify  two- 
phase  flow  patterns  [2].  In  addition,  encouraging  results  involving  the 
identification  of  two-phase  flow  patterns  of  refrigerants  have  been 
obtained  using  a  capacitance  probe  that  produced  void  fraction  signals 
[10].  However,  a  single  capacitance  probe  (as  in  [2,10])  may  not  fully 
characterize  the  liquid-vapor  spatial  distribution  within  a  flow  passage. 
This  could  lead  to  incorrect  assumptions  about  the  actual  flow  behavior. 
In  contrast,  electrical  capacitance  tomography  (ECT)  can  provide  a 
nearly  instantaneous  view  of  the  liquid-vapor  distribution  within  the 
system  without  optical  access.  Tomography  may  be  used  to  derive  the 
permittivity  distribution  from  capacitance  data  [8,9].  The  permittivity 
distribution,  in  turn,  can  provide  the  spatial  distribution  of  the  liquid 
and  vapor  phases.  Many  past  research  efforts  involving  ECT  have 
focused  on  industrial  applications  where  qualitative  results  were 
sufficient  [8].  The  use  of  ECT  in  detailed  studies  of  liquid  and  vapor  in 
horizontal  flow  has  been  limited,  much  less  with  the  use  of  dielectric 
refrigerants  which  are  of  interest  here  [8,10,1 1  ]. 

Artificial  neural  networks  are  used  for  pattern  recognition  and  trend 
prediction  involving  complex  processes.  In  an  artificial  neural  network, 
the  neurons  (often  called  nodes)  receive  input  signals,  and  each  node 
calculates  an  individual  output  using  a  weighted  sum  and  nonlinear 
activation  function.  Learning  is  achieved  by  the  adjustment  of  these 
weights  [12].  Past  studies  suggest  that  there  is  potential  for  using  a 
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Nomenclature 

Ai 

Individual  pixel  area  [m^] 

Aj 

Sum  of  all  pixel  areas  [m^] 

D 

Tube  diameter  [m] 

Db 

Bubble  diameter  [m] 

ECT 

Electrical  capacitance  tomography 

L 

Length  [m] 

PDP 

Probability  density  function 

t 

Time  [s] 

yM 

Center  of  mass  vertical  coordinate  normalized  by  the 
tube  diameter 

yc(t)pDF 

Probability  density  function  for  center  of  mass  (normal¬ 
ized)  vertical  coordinate 

yi 

Vertical  coordinate  normalized  by  the  tube  diameter  for 
the  i^^  pixel 

six,  t) 

Permittivity  determined  by  tomography  [E/m] 

Liquid  permittivity  [P/m] 

Vapor  permittivity  [P/m] 

e*(x, t) 

Normalized  permittivity 

e*(t) 

Spatial  average  of  normalized  permittivity 

®*(bpDf 

Probability  density  function  of  ^(t) 

e?(x,  t) 

Normalized  permittivity  for  an  individual  pixel 

mt)) 

Time  averaged  ^(t) 

{s*  (t)  Kurtosis  of  time  averaged  s*  (t) 

{^it))sKEw  Skewness  of  time  averaged  ^(t) 

Variance  of  time  averaged  s*  (t) 

X 

Spatial  locations,  x  and  y  [m] 

neural  network  to  objectively  classify  liquid-vapor  distribution  data 
[13-15].  These  studies  have  used  measured  or  simulated  impedance 
(conductance)  signals  rather  than  input  from  ECT  [13-16].  In  addition, 
nearly  all  previous  flow  identification  studies  that  have  used  neural 
networks  were  performed  for  vertical  flows. 

In  this  paper,  the  use  of  ECT  together  with  neural  networks  to 
identify  liquid-vapor  flow  patterns  is  explored.  Eor  this  purpose,  exper¬ 
iments  involving  the  horizontal  flow  of  refrigerant  R-134a  through  a 
tube  of  small  diameter  (7  mm)  were  conducted.  A  horizontal-two 
phase  flow  can  be  categorized  into  one  of  several  flow  patterns  which 
may  include  bubbly,  plug,  slug,  stratified-wavy,  and  annular  flows 
[17].  Here,  previous  work  is  extended  by  including  additional  transi¬ 
tional  flow  patterns  corresponding  to  four  of  the  above  flow  patterns 


for  horizontal  flow.  The  names  of  the  transitional  flow  patterns  are 
derived  from  these  four  flow  patterns  and  are  given  here  as  the 
bubbly-transitional,  plug-transitional,  slug-transitional,  and  stratified- 
wavy-transitional  patterns.  In  the  current  work,  high-speed  images  of 
the  flow  patterns  were  recorded  for  purposes  of  initial  human  classifica¬ 
tion  and  final  training  verification.  Processed  permittivity  data  obtained 
from  two-dimensional  tomograms  was  used  to  train  neural  networks.  A 
goal  was  to  classify  nine  horizontal  two-phase  flow  patterns  with 
reasonable  speed  using  neural  networks  as  a  predictive  tool.  Input 
information  for  the  neural  networks  included  the  spatially  averaged 
permittivity,  center  of  mass  location,  and  their  probability  density 
functions  (in  time).  It  also  included  four  statistical  moments  (in  time) 
for  spatially  averaged  permittivity  data.  Another  objective  was  to  deter¬ 
mine  which  subsets  of  input  data  provide  high  identification  success  for 
the  flow  patterns  when  analyzed  by  a  neural  network. 

2.  Experimental  setup 

To  explore  the  use  of  ECT  in  the  identification  of  two-phase  flow  pat¬ 
terns,  laboratory  experiments  were  performed  in  which  liquid-vapor 
flow  patterns  were  generated  for  flow  in  a  horizontal  tube.  Eig.  1 
shows  a  schematic  of  the  experimental  arrangement  in  which  liquid 
R-134a  was  pumped  through  a  heater  to  produce  two-phase  flow.  To 
obtain  different  liquid-vapor  flow  patterns,  the  volumetric  flow  rates 
were  adjusted  in  the  range  0.1  to  0.5  L/min,  while  varying  the  heater 
power  between  0  and  500  W.  Downstream  from  the  heated  section, 
the  two-phase  flow  entered  a  fused  quartz  observation  section  (tube 
with  7  mm  ID).  The  observation  section  permitted  imaging  using  a 
high  speed  video  camera  (Phantom  V4.2)  and  had  thermocouples  and 
pressure  transducers  located  at  1.2  m  increments  along  its  length.  The 
high  speed  camera  was  used  to  compare  actual  images  of  liquid-vapor 
flow  with  ECT  characterization.  R-134a  passed  through  the  cylindrical 
ECT  sensor  (ITS,  0690).  Eig.  1  shows  one  flow  path  for  R-134a  and 
another  for  cooling  water  which  was  used  to  condense  R-134a  vapor 
in  the  condenser  (Lytron,  LL510G02).  The  water  was  cooled  by  a  chiller 
(PolyScience,  4260  T). 

Table  1  lists  the  measurement  uncertainties  associated  with  the 
thermocouples,  pressure  transducers  (Omega  PX-409),  and  flow 
meter  (McMillan)  shown  in  Pig.  1.  Experiments  were  performed  with 
the  refrigerant  saturation  conditions  at  the  ambient  temperature  of 
20  °C.  Relatively  small  temperature  differences  (~1  °C)  between  the 
ambient  and  refrigerant  in  the  observation  section  of  tubing  provided 
reasonable  grounds  to  neglect  the  heat  transfer  between  them. 

The  ECT  system  consists  of  a  multi-electrode  sensor,  electronics  for 
capacitance  determination,  and  data  acquisition  components  [18].  The 
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Table  1 

Measurement  uncertainties. 


Instrument 

Type 

Calibration  range 

Uncertainty 

Thermocouple 

TypeT 

5-100 °C 

±0.1  °C 

Pressure  transducer 

Gauge 

0-1034  kPa 

±0.3  kPa 

Flow  meter 

Turbine 

0.10-1.00  L/min 

±  0.005  L/min 

sensor  fits  around  a  section  of  the  flow  passage  and  measures  changes 
in  capacitance  values  between  electrode  pairs.  Fig.  2a  shows  the  config¬ 
uration  for  eight  electrodes  (51  mm  length),  insulation,  grounded 
screen,  and  measurement  volume.  Starting  with  electrode  one,  a  voltage 
signal  was  applied  to  it.  The  capacitance  was  then  measured  from  the 
remaining  electrodes,  two  through  eight.  Next,  electrode  two  became 
electrified,  and  electrodes  three  through  eight  were  used  to  measure 
the  capacitance.  This  process  was  repeated  around  the  sensor  until  all 
(28)  pairs  of  measurements  were  completed.  A  single  tomogram  was 
generated  from  these  28  capacitance  measurements.  In  addition,  a 
sensing  unit  conditioned  the  sensor  signals  and  transmitted  them  to  a 
data  acquisition  system.  Details  of  the  hardware  and  experimental 
procedures  are  available  elsewhere  [19]. 

Each  tomogram  provided  a  two-dimensional  representation  of  the 
permittivity  obtained  from  capacitance  values  at  each  sample  time. 
Since  the  internal  sensor  volume  encompassed  a  portion  (51  mm)  of 
the  flow  path,  the  two-dimensional  permittivity  values  were  calculated 
from  a  three-dimensional  measurement  volume  and  averaged  in  time. 
Thus,  the  sensor  length  along  the  flow  path  as  well  as  the  sampling 
rate  were  practical  limits  for  the  sensor  resolution.  At  the  sampling 
rate  of  ~48  Hz,  20  ms  was  necessary  to  acquire  the  data  necessary  for 
one  two-dimensional  tomogram  with  a  delay  of  1  ms  between  data 
sets.  In  this  period,  a  two-phase  feature  with  a  speed  of  ~1  m/s  would 
be  detected  by  the  sensor  two  times  before  leaving  the  sensing  area. 
In  the  present  experiments,  the  greatest  velocity  of  a  vapor  feature 
was  -0.85  m/s.  Thus,  the  sampling  rate  was  reasonable. 

Eq.  ( 1 )  below  was  used  to  normalize  the  permittivity^ (x,  t),  that  was 
determined  by  tomographic  reconstruction.  Here,  x  represents  the  spa¬ 
tial  variables,  x  and  y,  and  t  is  time.  This  normalization  assigns  s*  (x,  t)  to 
be  zero  for  only  vapor  (Cg)  and  one  for  liquid  (Cf)  alone.  Since  non-linear 
electric  field  responses  may  generate  values  outside  the  normalized 


range,  a  baseline  reference  value  (Cf  or  8g)  was  recorded  while  the 
sensor  contained  only  liquid  (8f)  or  only  vapor  (8g).  A  tomogram  may 
have  regions  that  are  assigned  values  of  8*(x,  t)  between  zero  and  one 
since  a  given  pixel  may  represent  a  liquid-vapor  mixture. 

e*(x,t)  =  (e(x,t) -%)/(£/-%)  (1) 

In  this  work,  the  linear  back  projection  algorithm  was  used  in  the 
software  provided  by  the  instrument  manufacturer  for  the  tomographic 
reconstruction.  This  two-dimensional  (cross-sectional)  representation 
is  proportional  to  8*(x,  t)  within  the  sensor  volume  [18].  There  are  812 
pixels  which  comprise  the  computational  domain  for  the  reconstruction 
process  (Eig.  2b).  Since  the  cross-sectional  area  is  circular,  some  ele¬ 
ments  were  outside  the  circular  boundary  and  not  used.  In  addition, 
sensitivity  maps  generated  by  the  manufacturer  were  used  to  linearize 
the  measured  capacitance  values  which  were  ultimately  used  in  the  re¬ 
construction  algorithm  [18].  Subsequent  processing  of  the  tomograms 
was  performed  and  is  described  in  the  next  section  of  this  paper. 

Table  2  lists  the  flow  patterns,  heating  rates,  volumetric  flow  rates 
(at  the  pump  outlet),  pressure  range,  tube  diameter,  and  time  duration 
for  each  recorded  flow  pattern  case.  The  experiments  were  conducted 
by  holding  the  volumetric  flow  rate  constant  and  incrementally  increas¬ 
ing  the  heater  power  (Eig.  3).  The  heater  power  was  held  at  a  given 
value  until  the  liquid-vapor  behavior  was  steady  or  steady-periodic. 
Visual  observations  of  the  liquid-vapor  behavior  were  made  while 
recording  ECT  data.  A  total  of  546  heater  power  and  volumetric  flow 
rate  combinations  were  examined. 

3.  Data  processing  and  analysis 

Tomography  was  used  to  obtain  £*(x,  t)  from  the  capacitance  mea¬ 
surements  which,  in  turn,  provided  the  two-dimensional  distribution 
of  the  liquid  and  vapor  phases  within  the  section  of  tubing  encompassed 
by  the  sensor.  One  issue  concerning  the  tomographic  data  was  to  deter¬ 
mine  an  acceptable  number  of  input  data  to  train  the  neural  network  in 
a  reasonable  time  period  (less  than  10  min  on  a  Dell  Precision  T7500 
workstation,  64  bit  dual  processor).  The  required  computational  time 
was  prohibitively  large  for  directly  using  all  812  pixel  values  from 
each  tomogram  (500  tomograms  were  recorded  for  each  run)  for  each 


Fig.  2.  Drawings  of  the  a)  cross-sectional  view  of  the  sensor  with  eight  electrodes,  b)  tomography  computational  domain  consisting  of  812  pixels,  and  c)  instantaneous  normalized  vertical 
coordinate  (yc(t))  corresponding  to  the  refrigerant  center  of  mass.  The  normalized  coordinates  (Xc,  Yc)  designate  the  origin. 
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Table  2 

Summary  of  two-phase  flow  experiments. 


Common  flow  patterns 
Transitional  flow  patterns 

Heater  power  (W) 

Heat  flux,  q"  (W/cm^) 
Volumetric  flow  rate  (L/min) 
Pressure  (kPa) 

Tube  diameter  (mm) 
Duration  of  each  observed 
flow  pattern  case  (s) 


Bubble,  plug,  slug,  strategy  wavy,  annular 

Bubbly  transitional,  plug  transitional,  slug  transitional, 

stratified  wavy  transitional 

10-300 

0.40-12 

0.10-0.50 

500-800 

7 

30 


of  the  546  experimental  runs  as  input  data  (812  x  500  x  546  ~  2.2  x  10^ 
data  points)  to  train  a  neural  network.  It  was  proposed  that  a  combina¬ 
tion  of  spatial-  and  time-  averaging  together  with  other  quantities 
might  allow  for  rapid  training  of  accurate  neural  networks.  Here,  we 
explore  the  use  of  the  following  as  input  data  for  neural  network 
flow  pattern  identification:  the  spatial  average  of  £*(x,  t)  (represented 
by  ^(t)),  the  probability  density  function  of  ^(t),  the  refrigerant  local 
center  of  mass  (i.e.,  the  vertical  coordinate  normalized  by  the  tube 
diameter  and  denoted  by  3/c(t)),  and  the  probability  density  function 
of  yc(t).  In  addition,  the  time  average  of  ^(t)  and  the  associated 
variance,  kurtosis,  and  skewness  values  are  considered. 

^(t)  is  given  by  Eq.  (2),  where  Ai  is  the  individual  pixel  area  (Fig.  2b), 
£•  (x,  t)  is  the  normalized  permittivity  for  an  individual  pixel.  At  is  the 
sum  of  all  (812)  pixel  areas,  and  t  is  time. 

-|  812 

(2) 

Ai  in  Eq.  (2)  is  identical  for  all  pixels  and,  thus,  ^(t)  was  obtained  by 
taking  the  arithmetic  mean  of  the  s-  (x,  t)  values  for  each  tomogram  at  a 
particular  time,  t.  Since  £*(x,  t)  is  unity  for  the  presence  of  liquid  alone 
and  zero  for  vapor  alone,  ^(t)  represents  the  time-varying  liquid 
fraction  in  the  sensor  volume.  Consequently,  1  —  ^(t)  represents  the 
time-varying  void  fraction. 

For  each  flow  pattern  test  point  (i.e.,  each  heater  power  and 
volumetric  flow  rate  combination),  500  tomograms  were  generated 
over  a  sampling  period  of  30  s.  The  probability  density  function  (PDF) 
provided  a  way  to  observe  the  range  of  ^(t)  occurring  with  time  for  a 


given  case.  Here,  the  PDF  of  ^(t)  is  denoted  by  ^{t)pjjp.  To  estimate  ^ 
(t)pT,f ,  MATLAB  was  used  to  first  produce  a  histogram  by  separating 
the  500  values  of  ^(t)  (obtained  from  the  500  tomograms)  into  bins 
[20].  Preliminary  work  showed  that  calculated  histograms  were  sensi¬ 
tive  to  bin  sizes  below  25  but  were  insensitive  to  those  in  the  range 
25-100.  Thus,  50  bins  were  used  for  determining  all  histograms.  The 
estimate  of was  then  generated  by  a  smooth  fit  to  the  histogram. 
In  addition,  MATLAB  was  used  to  estimate  the  time  average,  {^(t)),  var¬ 
iance,  {^(t))v7y?.  skewness,  {^(t))s/c£w-  kurtosis,  {^{t))jajRT  of  ^(t) 
using  values  from  the  500  tomograms,  (^(t))  provides  a  representative 
value  for  the  time- varying  void  fraction,  {^(f))v7y?  provides  an  indication 
of  the  fluctuation  amplitude  of  {^(t)),  and  {^{t))sKEw  provides  a  mea¬ 
sure  of  the  asymmetry  of  the  probability  density  function.  In  terms  of 
a  physical  interpretation,  a  time  series  that  has  numerous  intermittent 
extreme  events  would  have  a  relatively  high  Thus,  a  high 

kurtosis  would,  for  example,  be  associated  with  the  rapidly  alternating 
occurrence  of  large  vapor  bubbles  and  liquid. 

A  running  average  was  used  to  estimate  the  number  of  data  samples 
required  to  obtain  a  stable  time  average  [21  ].  The  running  average  of 
^(t)  was  used  to  determine  an  acceptable  number  of  data  points  to 
estimate  {^(t)).  When  the  fluctuations  in  the  running  average  were 
within  a  +  5%  band  of  the  stable  average,  the  time  average  was  assumed 
to  be  stable  and  the  number  of  data  points  sufficient. 

Since  the  vertical  distribution  of  the  liquid  and  vapor  phases  in  a  hor¬ 
izontal  flow  can  be  influenced  by  gravity,  it  may  be  advantageous  to  use 
the  vertical  coordinate  of  the  center  of  mass,yc(t),  of  the  tomogram 
(Fig.  2b)  as  a  neural  network  input.  Since  the  origin  of  the  coordinate 
system  is  at  the  center  of  the  tube,  3/c(  t)is  zero  for  a  single  phase  flow 
within  the  tube.  Because  the  positive  direction  for  the  y  axis  is  up  with 
respect  to  the  gravity  vector,  a  two-phase  flow  pattern  will  result  in 
negative  3/c(t)  values  due  to  gravity  drawing  liquid  to  the  bottom  of 
the  tube.  £•  (x,  t)  and  the  densities  of  the  saturated  vapor  (pg)  and  liquid 
(p/)  were  used  to  estimate  the  fluid  mass  associated  with  a  given  pixel 
for  the  center  of  mass  calculation  as  given  by  Eq.  (3).  In  Eq.  (3),yi  repre¬ 
sents  a  vertical  pixel  coordinate  that  is  normalized  by  the  tube  diameter. 

812  /812 
yc(f)  =  E  [  (ps + f)  {pf  -pg)  )3'i]  /  E  [pg + 4  (X- 1)  (pf-pg)  ] 

i=l  /  i=l 

(3) 
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Fig.  3.  Visually  observed  flow  patterns  as  a  function  of  flow  rate  and  heater  power. 
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Following  the  calculation  procedure  with  ^(t),  PDFs  for3/c(t)  were 
generated  and  are  referred  to  as  3/c(f)pDF- 

In  this  work,  the  MATLAB  Neural  Network  package  was  used  as  the 
computational  tool  for  identifying  flow  patterns.  A  feedforward  neural 
network  was  selected  for  use  due  to  its  simplicity  and  ability  to  general¬ 
ize  to  new  situations  [20].  This  particular  feedforward  neural  network 
had  a  hidden  layer  with  nodes  that  used  a  sigmoid  transfer  function 
[20].  Training  of  a  neural  network  was  performed  to  identify  the  con¬ 
stants  that  comprise  the  weights  and  biases  of  each  neuron.  The  number 
of  hidden  nodes  was  varied  with  different  training  sessions  to  deter¬ 
mine  an  optimum  number  which  would  permit  rapid  calculations 
while  avoiding  overfitting.  In  this  work,  50  was  the  maximum  number 
of  hidden  nodes  used.  A  single  output  layer  was  generated  as  a  floating 
point  number  and  rounded  to  the  nearest  integer  value,  corresponding 
to  a  flow  pattern  prediction. 

The  feedforward  network  together  with  the  Levenberg-Marquardt 
training  algorithm  provide  rapid  calculation  speed  and  are  known  for 
successful  pattern  identification  [20].  Training  was  terminated  when 
the  generalization  no  longer  improved  which  was  indicated  by  an  in¬ 
crease  in  the  mean  square  error.  The  neural  network  was  retrained  by 
supplying  it  with  additional  randomized  subsets  of  the  original  input 
and  target  data.  This  was  repeated  until  the  neural  network  predicted 
the  observed  flow  patterns  with  an  acceptable  level  of  success,  defined 
here  to  be  the  correct  identification  of  flow  patterns  in  at  least  90%  of  the 
testing  data.  In  the  present  work,  the  neural  network  was  retrained  up 
to  twenty-five  times  to  reach  an  acceptable  level  of  success.  The  data 
used  for  testing  were  not  part  of  the  network  training  and,  thus,  provid¬ 
ed  an  independent  indication  of  neural  network  performance.  That  is, 
how  well  the  network  generalized  beyond  the  supplied  data  [20].  Data 
from  546  experiments  was  randomly  selected  such  that  30%  were 
used  for  the  validation  and  testing  of  a  neural  network,  and  70%  were 
used  for  training. 

The  neural  network  input  included  combinations  of  ^(t), 

(e^(t)),  {^(t))vAR,  {^(t))sKEw<(£'‘(t))KURT<  ydt).  and  ydOpDF  derived 
from  the  tomograms.  The  number  of  values  used  for  each  neural  network 
input  depended  on  the  type  of  data  itself  The  number  of  values  used  for 
^(t),  ^(OpDf’-Vc(t),  and3/c(t)pDFare  described  in  the  next  section.  The 
quantities  (?(t)),  (^(t))v-^,  (^(t))sK£w.  and  (^(t)>KURT  were  each  a 
single  value;  they  introduced  four  values  for  network  training  for  each 
data  set. 

4.  Results  and  discussion 

In  this  section,  the  flow  patterns  are  first  characterized  by  high  speed 
imaging.  This  is  followed  by  a  characterization  of  flow  patterns  by  ^(t), 
^(f)pDF’  yc(f),  and  yc(t)pDF-  The  performance  of  the  neural  network 
inputs  in  classifying  two-phase  flow  patterns  is  then  described. 

4.L  High  speed  imaging 

High  speed  imaging  was  used  for  the  initial  human  classification  of 
the  flow  patterns.  Table  3  shows  representative  images  (view  perpen¬ 
dicular  to  the  tube  axis)  and  cross-sectional  drawings  (view  along  the 
tube  axis)  for  each  flow  pattern.  It  is  useful  to  examine  the  images  of 
Table  3  together  with  the  plot  in  Fig.  3,  which  shows  the  dependence 
of  the  flow  pattern  on  both  the  input  heater  power  and  volumetric 
flow  rate.  Generally,  there  were  clear  visual  demarcations  among  the 
flow  patterns.  However,  overlapping  liquid-vapor  behavior  existed 
due  to  the  stochastic  nature  of  the  transition  patterns.  Care  was  taken 
to  label  the  transitions  correctly,  or,  in  the  case  of  great  uncertainty, 
eliminate  them  altogether.  Lastly,  Table  3  lists  the  time-averaged  void 
fraction  range  for  each  flow  pattern.  As  described  earlier,  the  void 
fraction  estimate  is  given  by  1  —  {^(t)). 

Table  3  shows  a  representative  image  of  bubbly  flow.  Bubbly  flow 
was  the  first  flow  pattern  encountered  after  single-phase  liquid  flow 


as  the  heater  power  was  increased  from  zero  with  the  volumetric  flow 
rate  at  the  highest  setting  (Table  2,  0.5  L/min).  As  illustrated  in  Table  3 
(drawing  of  view  along  tube  axis),  discrete  bubbles  with  diameters 
(Db)  smaller  than  the  tube  radius  (D/2)  moved  along  the  upper  tube 
surface  [17].  As  the  heater  power  was  further  increased  at  a  constant 
flow  rate  (Fig.  3),  bubbly  flow  evolved  into  what  is  termed  bubbly  tran¬ 
sitional  flow.  The  image  of  bubbly  transitional  flow  in  Table  3  is  similar 
to  that  of  bubbly  flow  but  elongated  bubbles  move  along  the  upper  tube 
surface.  The  drawing  along  the  tube  axis  indicates  that  Db  in  the  bubbly- 
transitional  pattern  were  smaller  than  D/2,  but  the  bubble  length 
(L)  could  approach  2D. 

Plug  flow  occurred  when  small  adjacent  bubbles  coalesced  into  larg¬ 
er  bubbles  separated  by  liquid  regions  [17].  In  addition,  Db  continued  to 
be  smaller  than  D,  but  L  was  several  D,  and  the  bubble  frequency 
increased  significantly.  Bubbles  in  plug-transitional  flow  (Table  3) 
have  L  and  Db  similar  to  those  in  plug  flow.  However,  plug-transitional 
bubbles  have  an  elongated  front  section  followed  by  a  tail  of  smaller 
diameter.  Slug  flow  evolved  from  plug  transitional  flow  where  L 
approached  SOD.  The  bubbles  had  a  consistent  bullet  shape  with  a 
rounded  nose  and  tail.  Also,  Db  began  to  approach  D.  Liquid  regions  sep¬ 
arated  the  vapor  bubbles  and  often  entrained  small  dispersed  bubbles 
[17].  Slug-transitional  flow  is  similar  to  slug  flow  in  that  the  bubbles 
have  a  bullet  shape  with  a  short  tail  (not  in  image).  However,  L  was 
greater  than  that  for  slug  flow  with  smaller  liquid  regions. 

Stratified  wavy  flow  occurred  at  high  heater  powers  (Fig.  3).  The 
image  of  stratified  wavy  flow  shows  that  liquid  slugs  no  longer 
occurred,  and  a  vapor  layer  existed  along  the  upper  portion  of  the 
tube.  Kelvin-Helmholtz  instabilities  are  responsible  for  waves  which 
propagated  along  the  tube  length,  often  contacting  the  upper  tube 
surface  [17,22].  As  in  stratified  wavy  flow,  interfacial  waves  were  also 
present  in  stratified  wavy  transitional  flow.  The  waves  were  more  fre¬ 
quent  but  their  crests  did  not  contact  the  upper  tube  surface.  However, 
they  rose  up  on  the  tube  wall  wetting  the  lateral  surface  [17,22].  Table  3 
shows  an  example  of  annular  flow  consisting  of  liquid  flow  enveloping  a 
flowing  vapor  core  with  small  waves  along  the  liquid-vapor  interface. 
Annular  flow  is  believed  to  occur  when  the  vapor  velocity  increases 
beyond  that  in  stratified  wavy  transitional  flow  [17,22].  The  liquid 
layer  is  thicker  at  the  bottom  due  to  gravity  and  did  not  always  extend 
all  the  way  to  the  top. 

4.2.  ^(t)  and  ^(t)pDp  characteristics  for  dijferent  flow  patterns 

The  flow  patterns  are  now  described  in  terms  of  ^(t)  and  ^{t)pjjp. 
Fig.  4  shows  representative  plots  of  ^(t)  with  time  together  with 
£*(f)pDF  for  nine  flow  patterns  in  the  order  of  increasing  void  fraction. 
The  plots  correspond  to  the  images,  flow  rates  and  heater  powers  of 
Table  3.  Bubbly  flow  (Fig.  4a,  flow  rate  of  0.5  L/min)  has  values  of  ^(t) 
in  the  range  0.95-1.  ^{t)pp,p  for  bubbly  flow  (Fig.  4b)  has  a  single  peak 
near  the  upper  end  of  the  permittivity  range,  indicating  that  most  of 
the  flow  is  liquid.  As  the  heater  power  increased  at  the  flow  rate  of 
0.5  L/min,  bubbly  transitional  flow  (Fig.  4c  and  d)  resulted  which  had  a 
broader  ^(t)poF  with  lower  values  of  ^(t).  At  a  flow  rate  of  0.3  L/min 
and  heater  power  near  30  W  (Table  3),  the  ^(t)  associated  with  plug 
flow  (Fig.  4e)  have  a  nearly  periodic  behavior  with  values  ranging  from 
0.56  to  0.98.  This  nearly  periodic  behavior  is  confirmed  by  the  image 
in  Table  3  of  flowing  liquid  plugs  that  are  regularly  separated 
by  vapor.  ^(t)poF  for  plug  flow  (Fig.  4f)  has  two  distinguishable  peaks 
for  ^(t)  of  0.78  and  0.93  due  to  the  presence  of  large  volumes  of  alter¬ 
nating  liquid  and  vapor  phases.  The  liquid  plugs  appear  less  periodic 
in  the  plug  transitional  flow  pattern  (Fig.  4g),  and  the  corresponding 
^{t)pDP  (Fig.  4h)  is  similar  although  broader  than  ^{t)pp,p  for  plug  flow. 
Fig.  4i  shows  that^(t)  of  the  slug  flow  pattern  has  large,  relatively  rapid 
fluctuations  that  are  associated  with  the  passage  of  long  vapor  bubbles 
and  short  liquid  slugs  (Table  3).  Fig.  4j  shows  that  the  ^{t)ppfp  for  slug 
flow,  has  a  broad  range.  For  slug  transitional  flow  (Fig.  4k),  the  heights 
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Table  3 

Representative  flow  pattern  images  and  characteristics. 
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of  the  peaks  in  ^(t)  have  generally  increased  and  the  ^(t)  values  in  the 
^(t)pDF  plot  range  from  0.21  to  0.98  (Fig.  41).  Relative  to  slug  flow,  this 
indicates  sequential  vapor  bubbles  with  larger  volumes  that  are  sepa¬ 
rated  by  a  greater  distances.  The  ^(t)  associated  with  stratified  wavy 
fiow  (Fig.  4m)  has  large  and  small  intermittent  peaks  relative  to  ^(t) 
of  the  slug  and  slug  transitional  patterns.  The  intermittent  peaks  are 
due  to  interfacial  waves  which  propagated  along  the  tube  length 
(Table  3).  Relative  to  stratified  wavy  fiow,  Fig.  4o  shows  that  ^(t)  for 
stratified  wavy  transitional  fiow  has  significantly  smaller  peaks.  In  addi¬ 
tion,  Fig.  4p  shows  that  ^(t)p£)p  is  entirely  below  0.5  which  also  means 
that  the  void  fraction  was  always  greater  than  0.5.  The  amplitude  of  the 
fiuctuations  in^(t)  for  annular  fiow  is  small  (Fig.  4q)  relative  to  those  of 
stratified  wavy  or  stratified  wavy  transitional  flows.  This  is  confirmed 
by  the  observation  that  interfacial  waves  were  no  longer  visible 


(Table  3).  In  addition,  the  plot  of  ^(t)pp,p  for  annular  fiow  has  a  peak 
with  no  tail  which  accompanies  a  large  void  fraction  (Table  3).  This 
might  be  anticipated  from  the  images  of  a  fiowing  vapor  core  and  liquid 
annulus  within  the  tube  (Table  3)  for  these  fiow  patterns. 

4.3.  yc(t) ,  and  yd  t)  pdf  Characteristics  for  different  flow  patterns 

For  a  liquid-vapor  fiow  that  is  horizontal  with  respect  to  gravity, 
the  gas  phase  is  near  the  upper  surface  of  the  fiow  passage  with  the 
liquid  phase  below  if  the  tube  is  long  enough  for  the  vapor  to  stratify. 
Consequently,  it  is  helpful  to  use  the  vertical  coordinate  of  the  tomo¬ 
gram  centroid  location  as  a  neural  network  input.  Fig.  5  shows 
yc(t),  and  yc(t)pDF  for  each  fiow  pattern.  For  reference,  a  y-coordinate 
value  of  0  corresponds  to  the  center  of  the  tube.  The  peaks  and  troughs 
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of  yc(t)  in  Fig.  5  generally  track  with  those  of  ^(t)  in  Fig.  4  for  the  corre¬ 
sponding  patterns.  This  occurs  for  two  reasons.  The  first  is  that  the 
center  of  mass  shifts  downward  (upward)  as  the  vapor  fraction 
increases  (decreases)  due  to  the  difference  in  phase  densities  in  a 


gravity  field.  The  second  reason  is  that  ^(t)  decreases  with  increasing 
vapor  volume.  Figs.  4  and  5  show  that^(t)p^p  and  yc(t)pDF  both  broaden 
beyond  a  single  peak  and  then  ultimately  return  to  a  single  peak  from 
bubbly  to  annular  flows. 
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Fig.  5.  yc(t)  and  yc(t)pDF  behavior  for  bubbly  (a  and  b),  bubbly  transitional  (c  and  d),  plug  (e  and  f),  plug  transitional  (g  and  h),  slug  (i  and  j),  slug  transitional  (k  and  1),  stratified  wavy 
(m  and  n),  stratified  wavy  transitional  (o  and  p),  and  annular  (q  and  r)  flows. 


4.4.  Neural  network  studies 

It  is  useful  to  determine  which  combinations  of  input  data  types 
derived  from  ECT  tomograms  allow  rapid,  classification  of  two-phase 
flow  patterns  by  a  feedforward  neural  network.  The  current  study  was 


carried  out  in  two  stages  with  increasing  complexity.  In  the  first,  bubbly, 
plug,  slug,  stratified-wavy,  and  annular  patterns  were  considered  since 
they  are  relatively  well  defined  and  have  been  treated  elsewhere  [17, 
22].  It  was  proposed  that  combinations  of  input  data  types  that  resulted 
in  high  success  rates  for  classifying  fiow  patterns  in  the  first  stage  could 
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then  be  used  in  the  second,  more  challenging  stage  to  identify  nine  flow 
patterns. 

Table  4  compares  the  classification  success  rates  for  eight  neural 
network  input  vectors  for  the  five  basic  flow  patterns  used  in  the  first 
stage  of  this  study.  The  first  column  provides  the  flow  pattern  category, 
and  the  second  lists  the  number  of  experimental  runs  for  each  flow 
pattern.  The  third  through  the  tenth  columns  provide  the  classification 
success  rates  for  input  that  was  comprised  of  individual  and  combined 
data  types.  Here,  the  success  rate  is  defined  as  the  percentage  of  flow 
patterns  predicted  correctly.  Based  on  past  work  involving  vertical 
flows  and  other  instrumentation,  there  was  reason  to  believe  that  the 
use  of  neural  network  input  values  derived  from  time  averaging  or 
probability  density  functions  have  potential  for  reasonable  success 
rates  [13,14].  Thus,  the  seventh  column  lists  the  success  rates  for  flow 
pattern  identification  that  used  {^(t)),  {^{t))vAR^  (^i^))sKEw^ 
{^i^))i<uRH  together  as  input  for  each  experiment.  Similarly,  columns 
8  through  10  list  the  success  rates  for  input  comprised  of  (^(t)), 
(  ^*(0)sk£W5  {^*W)kurt  with  £*(t)pj3p  and/or  y(.(t)pp,p  for 
each  experimental  run.  The  second  row  in  Table  4  provides  the  number 
of  input  data  points  for  each  data  type.  The  bottom  row  lists  the  average 
success  rate  of  the  neural  network  resulting  from  the  use  of  different 
input  data  types. 

Table  4  shows  that  using  either^(t)  or  yc(t)  alone  as  an  input  result¬ 
ed  in  the  lowest  overall  success  rates  with  the  least  success  in  defining 
slug  flow.  The  use  of  other  input  types  provided  success  rates  greater 
than  90%  for  all  flow  patterns.  Since  a  success  rate  of  90%  or  greater 
is  a  goal.  Table  4  shows  that  it  is  acceptable  to  use  any  listed  input 
other  than  ^(t)  or  yc(t)  to  identify  the  basic  flow  patterns.  The  average 
classification  success  rate  over  all  patterns  for  columns  5  and  7  shows 
that  the  first  four  statistical  (time)  moments  involving  ^(t)  provide 
information  that  is  only  slightly  better  than  that  of^(t)pp,p.  This  supports 
the  conclusion  of  Rosa  et  al.  [13]  that  the  first  four  statistical  (time) 
moments  involving  ^(t)  provide  information  essentially  equivalent  to 
that  of^(t)pp,p. 

The  bubbly,  plug,  slug,  and  stratified-wavy  flow  patterns  were  then 
subdivided  into  categories  which  included  transition  patterns.  In  these 
transition  regions  the  flow  may  oscillate  between  the  adjacent  patterns 
making  classification  difficult.  Consequently,  Table  5  shows  neural  net¬ 
work  classification  results  for  nine  flow  patterns.  The  number  of  runs 
(second  column)  listed  under  the  original  and  transition  flow  patterns 
in  Table  5  sum  to  the  original  number  of  runs  for  the  basic  flow  pattern 
in  Table  4  (i.e.,  the  transitional  flow  patterns  in  Table  5  are  subsets  of  the 
basic  flow  patterns  in  Table  4).  Lastly,  the  neural  network  input  types  of 
Table  5  are  the  same  as  those  in  Table  4  (i.e.,  columns  3-10). 


Table  5  shows  that  the  lowest  success  rates  occurred  using  either 
^(t)  or  yc(t)  alone  as  was  shown  in  Table  4.  In  fact,  the  success  rate 
using  either  as  a  network  input  was  zero  for  bubbly  and  bubbly  tran¬ 
sitional  flows.  Zero  success  rates  for^(t)  and  yc(t)  are  believed  to  be 
due  to  the  challenge  of  differentiating  among  the  strong  similarities 
between  bubbly  and  bubbly  transitional  patterns  (Figs.  4a,  c,  5a,  and 
c).  Table  5  shows  that  the  use  of  ^(t)p^p,  yc(t)pDF,  or  the  group  {^(t)) , 
{^W)vAR^  {^(0)s/c£w.  ^nd  {^(t))/ct/i?T’  input  had  similar  average 
overall  classification  success  rates  (last  row)  and  offered  significant  im¬ 
provement  in  classification  success  over  that  of  ^(t)  or  yc(t)  for  individ¬ 
ual  flow  patterns.  This  can  be  attributed  to  more  distinct  differences 
among  the  PDFs  and  time  averaged  statistics  for  different  patterns 
(PDFs  in  Figs.  4  and  5  and,  for  example,  the  time  averaged  statistics  as 
represented  by  the  void  fraction  in  Tables  3, 1-  ^(t).  The  bottom  row 
of  Table  5  shows  that  the  average  classification  success  rate  over  all 
flow  patterns  increases  when  the  time  averaged  statistics  input,  {^(t)), 
{^(f))sK£w.  ^nd  {^(t))j^pp,  is  combined  with  either  of  the 
input  values,  ^(f)poF  or  yc(t)pDF-  Further,  when  ^(t)poF  and  yc(t)pDF  are 
combined  with  time  average  statistics  into  a  single  input  (column  10), 
classification  success  rates  greater  than  90%  were  obtained  for  all 
patterns.  In  addition,  the  average  classification  success  rate  over  all 
flow  patterns  was  98.1%. 

The  classification  success  rates  for  the  above  nine  flow  patterns 
obtained  using  the  combination  of  time  averaged  statistics,  and 

yc(t)pDF  (i.e.,  column  10)  compare  well,  if  not  better,  than  those  of  past 
studies  involving  electrical  capacitance  measurements  (ECT  sensors) 
and  horizontal  flows.  For  example,  Jeanmeure  et  al.  [23]  used  an  identifi¬ 
cation  tree  to  classify  horizontal  flows  but  only  used  three  patterns.  Later, 
Yan  et  al.  [24]  used  parameters  obtained  from  steady,  two-dimensional 
numerical  simulations  of  nine  flow  patterns  to  train  a  back-propagation 
neural  network.  Their  numerical  simulations  were  correlated  with  exper¬ 
iments  involving  tomography  but  did  not  include  the  inherent  unsteady 
nature  of  an  actual  two-phase  flow.  In  addition,  two  of  their  patterns,  a 
tube  filled  with  liquid  or  empty,  were  relatively  simplistic. 

Table  5  lists  the  classification  success  rates  for  the  flow  patterns 
but  does  not  specifically  characterize  the  classification  errors  for  the 
combined  input  (e^(t)),  (^{t)}sKEw^  (^(t)}KURT-  ^(Opof 

and  yc(t)pDF-  An  error  value  within  ±0.5  of  zero  represents  correct 
identification,  otherwise  a  value  outside  these  bounds  is  considered  to 
be  incorrect.  For  a  total  of  546  experiments,  a  small  number  (14)  were 
incorrectly  classified  by  one  flow  pattern  and  a  much  smaller  number 
( 1 )  were  incorrectly  classified  by  two  patterns.  More  specifically,  99% 
of  the  experimental  runs  were  either  correctly  classified  or  misclassified 
by  no  more  than  one  flow  pattern.  Tables  4  and  5  together  show  that  the 


Table  4 

Comparison  of  Classification  Success  Rates  Using  Eight  Neural  Network  Input  Vectors  for  Five  Basic  Flow  Patterns. 


Flow  Pattern 

Number  of 

Runs 

s*(t) 

yc(t) 

£*(t)pDF 

Ycft)  PDF 

<6*(t)> 

<e*(t)> 

(e*(t)> 

(fi  (t))  VAR 

(t))  VAR 

(t))  VAR 

(t))  VAR 

(^*(t)}  SKEW 

SKEW 

(^*(t))  SKEW 

SKEW 

(fi  (t)}iajRT 

(fi  (t))KURT 

(X))  lOJRT 

i^))  lOJRT 

Ycft)  PDF 

£  (t)  PDF 

£  (t)  PDF 

Vcft)  PDF 

Number  of  Data  in  Input  Vector  for  Each  Run 

250 

250 

50 

50 

4 

54 

54 

104 

Bubbly 

102 

100 

100 

100 

100 

100 

100 

100 

100 

Plug 

98 

94.9 

96.9 

100 

100 

100 

100 

100 

100 

Slug 

129 

86.1 

86.1 

98.4 

99.2 

100 

99 

100 

100 

St.  Wavy 

129 

94.6 

89.1 

100 

100 

99.2 

100 

100 

100 

Annular 

88 

88.6 

95.5 

97.7 

100 

98.9 

100 

100 

100 

Average  Classification  Success  Rate  over  All  Patterns 

92.8% 

93.5% 

99.2% 

99.8% 

99.6% 

99.8% 

100% 

100% 
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Table  5 

Comparison  of  Classification  Success  Rates  Using  Eight  Neural  Network  Input  Vectors  for  Nine  Flow  Patterns. 


Flow  Pattern 

Number  of 

Runs 

e*{t) 

yc(t) 

£  (t)pDF 

yc(t)  PDF 

{s*it)) 

{s*it)) 

{s*it)) 

(e*(t)> 

(t))  VAR 

(t))  VAR 

(fi  (t))  VAR 

(fi  (t))  VAR 

(£*(t))sK£W 

SKEW 

(fi*(t))sKEW 

(fi*(t))sKEW 

(X))  lOJRT 

(fi  it))  KURT 

(fi  it))  KURT 

1 

1  * 

yc(t)  PDF 

PDF 

£*(t)  PDF 
yc(t)  PDF 

Number  of  Data  in  Input  Vector  for  Each  Run 

250 

250 

50 

50 

4 

54 

54 

104 

Bubbly 

52 

0 

0 

94.2 

94 

94.3 

100 

98.1 

100 

Bubbly  -  T 

50 

0 

0 

88 

96 

76.1 

98 

96.0 

94.0 

Plug 

58 

32.8 

31.1 

100 

97 

100 

100 

98.3 

100 

Plug  -  T 

40 

20.0 

22.5 

95 

98 

90.1 

100 

100 

100 

Slug 

89 

29.2 

41.6 

80 

79 

94.4 

93.3 

92.3 

96.6 

Slug-T 

40 

17.5 

35 

65 

63 

57.5 

85 

92.5 

100 

St.  Wavy 

69 

30.4 

24.6 

52.2 

71 

81.2 

84.1 

89.4 

97.0 

St.  Wavy-  T 

60 

13.3 

5 

67 

63 

81.7 

91.7 

90.0 

95.0 

Annular 

88 

39.8 

28.4 

95.5 

86 

98.9 

97.7 

98.9 

100 

Average  Classification  Success  Rate  over  All  Patterns 

20.3% 

20.9% 

81.9% 

82.9% 

86.0% 

94.4% 

95.1% 

98.1% 

In  this  table,  T  in  the  pattern  name  represents  the  word  “transitional.” 


input  combination  of  time  averaged  statistics,  £*{t)pjjp,  and  yc(t)pDF  can 
be  used  to  reasonably  classify  flow  patterns  observed  in  horizontal  flow. 

5.  Conclusions 

This  work  investigated  the  use  of  ECT  and  neural  networks  to 
identify  two-phase  flow  patterns  for  refrigerant  R-134a  flowing  in  a 
horizontal  tube.  Laboratory  experiments  were  performed  to  generate 
different  liquid-vapor  flow  patterns.  High-speed  images  were  recorded 
for  human  classification  of  the  flow  patterns.  The  corresponding  permit¬ 
tivity  data  obtained  from  tomograms  was  used  to  train  feedforward 
neural  networks  to  recognize  flow  patterns.  A  goal  was  to  classify  nine 
horizontal  two-phase  flow  patterns  using  neural  networks  as  a  predictive 
tool.  Another  was  to  determine  which  subsets  of  input  data  provide  high 
identification  success  for  the  flow  patterns  when  analyzed  by  a  neural 
network.  Previous  work  involving  horizontal  flows  was  extended  by 
considering  transitional  flow  patterns  associated  with  four  common 
flow  patterns  (within  a  small-diameter  tube).  A  unique  feature  of  the  cur¬ 
rent  work  is  the  use  of  the  vertical  center  of  mass  coordinate  in  pattern 
classification  for  horizontal  flow. 

The  highest  classification  success  rate  occurred  with  the  use  of  neu¬ 
ral  network  input  consisting  of  yc(t)pDF,  (^(0). 

(^{t))sKEw^  ^(0/<uRT*  The  combination  of  these  data  resulted  in  a 

1 00%  success  rate  for  the  basic  five  flow  patterns  and  a  high  average  suc¬ 
cess  rate  of  98.1%  for  nine  flow  patterns,  which  include  the  transitional 
patterns.  In  addition,  99%  of  the  experimental  runs  were  either  correctly 
classified  or  misclassified  by  only  one  flow  pattern.  Based  on  the  success 
rate  of  flow  pattern  identification  with  relatively  small  network  training 
and  run  times,  there  is  potential  for  the  use  of  ECT  combined  with  neural 
network  analysis  to  automatically  identify  two-phase  patterns  for 
horizontal  flow  in  thermal  process  control  applications. 
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